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Abstract 

A quantum gate is realized by specific unitary transformations operating on states representing 
qubits. Considering a quantum system employed as an element in a quantum computing scheme, 
the task is therefore to enforce the pre-specified unitary transformation. This task is carried out 
by an external time dependent field. Optimal control theory has been suggested as a method 
to compute the external field which alters the evolution of the system such that it performs 
the desire unitary transformation. This study compares two recent implementations of optimal 
control theory to find the field that induces a quantum gate. The first approach is based on the 
equation of motion of the unitary transformation. The second approach generalizes the state to 
state formulation of optimal control theory. This work highlight the formal relation between the 
two approaches. 

PACS number(s): 82.53.Kp 03.67.Lx 33.90.+h 32.80.Qk 
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I. INTRODUCTION. 



Quantum computation is based on implementing selected unitary transformations rep- 
resenting algorithms [|l[]. In many physical implementations the unitary transformation is 
generated using an external driving field. This driving field has to perform the quantum 
gate between, for example, two qubits, without altering the other levels which represent 
additional qubits. This means that the specific unitary transformation has to address a set 
of levels in an environment where other energy levels are present. An approach to correct 
this undesired coupling to other levels has been suggested for specific cases [0 but a general 
solution is not known. 

The presence of a large number of levels coupled to the external driving field is specially 
relevant in the implementation of quantum computing in molecular systems P, |], ||. Tesh 
et al. H proposed the use of optimal control theory (OCT) as a possible remedy. OCT 
0] is a well developed approach that allows to obtain the driving field which induces a 
desired transition between preselected initial and final states. However, the dependence on 
the particular transitions (characterized by the initial and final states) makes the traditional 
formulation of OCT inappropriate for that purpose. For example, if the unitary transfor- 
mation relates the initial states tpa- with the final states (the index k denotes all the 
relevant states involved in the transformation), the traditional OCT approach derives an 
optimal field e& for each pair {(fiik,ffk}- But the fields e& obtained are in general different, 
and then the transition induced by for the initial state (fin-' won't necessarily give the right 
final state <Pfy- To implement a given unitary transformation a single field e that relates all 
the relevant pairs {y^fc, is needed. 

An approach to overcome this problem was suggested in [0, |J. The idea is to generalize 
OCT to deal directly with the evolution operator. Recently, a different approach using 
simultaneous optimization of several state to state transitions has been suggested by Tesch 
and Vivie-Riedle ||. The purpose of this work is to compare the two approaches and two 
point the similarities and differences between them. This comparison has led to new insight 
into the use of optimal control theory for quantum systems. 

Unitary transformation optimization can be formulated in the following way. We consider 
a quantum system with a Hilbert space of dimension M, expanded by an orthonormal 
basis of states {\k)} {k = 1,...,M). In the following we will assume that the k states 
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correspond to the free Hamiltonian proper states. The selected unitary transformation 
is imposed on the subspace of the first iV energy levels of the system (N < M). For 
example, the N levels could correspond to the physical implementation of the qubit(s) 
embedded in a larger system. The additional levels (k = N + 1,...,M) are not directly 
involved in the target unitary transformation and they are generally considered as "spurious 
levels" coupled to the field. However, it is not always the case: an example is the proposal 
of implementing quantum computation using rovibronic molecular levels. In the simplest 
description two electronic surfaces are considered. Two rovibronic states of one of the 
electronic surfaces are chosen as the implementation of a qubit (in this case N = 2) and the 
unitary transformation is implemented using field induced transitions between the electronic 
surfaces. Out of the relevant subspace (k = 3,...,M) there are levels residing on both 
electronic surfaces, spurious in the sense that any leakage to them would destroy the desired 
final result, but at intermediate times they are used as intermediate storage which allow to 



carry out the desired unitary transformation || [HJ . 

The optimization objective is to implement a selected unitary transformation in the 
relevant subspace at a final time T. The target unitary transformation is represented by 
an operator in the system Hilbert space and denoted by O. For iV < M, the matrix 
representation of O in the basis {\k)} has two blocks of dimension N x N and (M — 
N) x (M — N). The elements connecting these blocks are zero. This structure means that 
population at the target time is not transferred between the two subspaces. Only the N x N 
block is relevant for the optimization procedure, and the other block remains arbitrary. 

The dynamics of the system is generated by the Hamiltonian H, 

H(t;e) = H -£e(t), (1) 

where Ho is the free Hamiltonian, e(t) is the driving field and ft, is a system operator 
describing the coupling (transition dipole operator in molecular systems). Eq. (|l|) can be 
generalized to more than one independent driving field, for example, controlling separately 
the two components of the polarization of an electro- magnetic field ||1 1|| . The dynamics of 
the system at time t is fully specified by the evolution operator U(t, 0; e). An optimal field 
e op t induces the target unitary transformation O on the system, at time T if 

U(T,0;e opi ) = e"^ 6 . (2) 
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Eq. (H) implies a condition only on the N x N block of the matrix representation of U. The 
phase 4>(T) is introduced to point out that in some cases the target unitary transformation 
O can be implemented up to an arbitrary global phase. The phase can be decomposed 
into two terms, 0i (T) + 2 (T) . The term (pi originates from the arbitrary choice of the origin 
of the energy levels which formally means that a term proportional to the identity operator 
can always we added to the Hamiltonian. The phase 0i is given by, 

Wr) = l^l (3) 

where Ek is the energy of the level k. By 02 we denoted other contributions to the global 
phase due to the structure of the unitary transformation and its arbitrariness for the levels 
k = N + 1,...,M. 



II. EVOLUTION EQUATION APPROACH. 



The optimization approach proposed in 
parameter r as 



is based on Eq. (0) by defining a complex 



N 

r(6;T; e ) = ^(A;|6 t U(T,0; e )|A;). 



(4) 



fe=i 



As O is a unitary transformation in the relevant subspace, r is a complex number inside a 
circle of radius N. Its modulus is equal to N only when the unitary transformation generated 
by the field U is equal to the target unitary transformation in the relevant subspace, except 
for a possible global phase. The modulus of r is a measurement of the fidelity of the target 
unitary transformation implementation by the field M. For iV = M, the sum in Eq. (ffl) is 
the trace of the operator product. 

In 0, U the optimization of the real part of r, or the imaginary part, or a linear com- 
bination of both was suggested as a method to find the optimal field. For simplicity we 
will consider the real case. The maximization of Re[r] can be formulated as a functional 
optimization. In this work we use the following form 

■ n I r T 

V7A;|6 T U(T,0;eo + Ae) |Jfe) -X |Ae| 

Jo 



J(U,B,Ae) = Re 



Re 



fe=i 

N 

E 

h=i 



2 dt 



(k\B[ ^ + lH(t ;eo + A e ) )V\k)dt 



(5) 



where U, B, eo, and Ae depends on t. The first term in the right-hand-side is the original 
objective. In our formulation, the other two terms are constrains depending on a reference 
field eo and a field correction Ae. The term including |Ae| 2 minimizes the total energy of 
the correction. The last term introduces the dynamics of the system under the field eo + Ae. 
B is an operator Lagrange multiplier and A a scalar Lagrange multiplier ||. The common 
OCT form for the functional || is recuperated setting eo = and interpreting Ae as the 
driving field. However, Eq. (|5|) offers some advantages in the interpretation of the equations 
derived from the functional. More elaborated constrains are possible [ F2| , |T3 |, but they are 
not relevant for the following discussion. 

Applying the calculus of variations, 5 J = 0, with respect to B, the Schrodinger equation 
for the evolution operator of the system is obtained, 

^ = -^H(t;e + Ae)U, (6) 

with the condition U(0, 0; eo + Ae) = 1. The variation of U give the Schrodinger evolution 
equation for the operator B, 

^ = -^H(t;6 + Ae)B t , (7) 

with the condition B^(T, T; eo + Ae) = O. B^ can be interpreted as the backwards propa- 
gation in time of the target unitary transformation O. It is related to U by 

B t (t,T;e) = U(t,T;e)6. (8) 

Eq. (§) and (0) represent the propagation forward and backwards in time of the boundary 
conditions of the problems, that is, the identity, 1, at time t = and the target unitary 
transformation, O, at time t = T. The variation of Ae leads to an equation for the correction 
to the field, 

1 N 

Ae(t) = -— - bn[ 5>l B(t, T; e + Ae) £ U(t, 0; e + Ae) \k)} . (9) 

k=l 

Eq. @, more than Eq. (|j), is the central result of the optimal control procedure and 
constitutes the starting point of the iterative algorithms devoted to determine the optimal 
field When an optimal field e opt is found Ae = 0, and from Eq. (|j) and Eq. (|9]), 

N 

Im[J2W&v\t,T;e opt )fiij(t,0-,e opt )\k)} =0. (10) 

k=l 
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Eq. (|T0| ) constitutes a condition for the optimal fields for the first approach. It is the 
base for the following analysis. Let us denoted by e a field that generated the target unitary 
transformation up to a global phase, U(T, 0; e) = e - ^ O. Using 

U(t,0,e) = U(t,r,e)U(T,0,e), (11) 

the diagonal block structure of the matrix representation of O in the basis {|/c)}, and the 
relation (k\p,\k) = 0, it is found that the left-hand-side of Eq. (0) gives 

N 

Im[ e ^( k \ U(*, T; e) 6 6 t tj\t, T; e) p, \k)} = . (12) 

k=l 

This result implies that any field inducing the target unitary transformation up to a global 
phase is a possible optimal solution of the optimization algorithm based on Eq. flj^). The 
convergence to such a solution will depend on the particular numerical implementation. 
It must be remarked that Eq. (pTOj) is a necessary but not a sufficient condition for the 
optimal field. For example, let us consider a target unitary transformation diagonal in the 
basis {\k)} denoted by Op. The unitary transformation generated by the free Hamiltonian 
XJ(ti,t2] e = 0) is also diagonal in that basis, and then, 

TV 

Im[^(A;| 6^U t (t,r;e = 0) £U(t,0;e = 0) \k)] =0. (13) 
fc=i 

However, in general U(T, 0; e = 0) is not equal to the target Od- These spurious solutions to 
the optimization can be avoided with a different choice of the initial guess for the algorithm. 

In optimization procedures based on Eq. ([]) the full operator propagation in Eq. (^) and 
(0) is not needed since only the action of U and O on the states \k) in the relevant subspace 
(k = 1, ...,N) appears. Then only the first N rows of the matrix operator representations 
are propagated. Denoting by U k ((B Jf ) h ) the k row of the matrix representation of U (B ) in 
the basis {|&)}, and being U* ((B^)*) the j element of the row (j = 1, M), the evolution 
equations @ and (0) take the form, 

dU k (t) 



at h 

d{B^f{t) i 



-H(t) U k (t) , 
H(t)(B^ k (t), (14) 



dt h 

with the conditions Uj(t = 0) = 5 jk and (B^) k (t = T) = O k respectively, being O k the 
matrix elements of the target unitary transformation O and H(t) the matrix representation 



of H. The 2N evolution equations (k = 1,...,N) are equivalent to the propagation of 
2N states of the system. The advantage is that in this case only 2(iV x M) elements are 
propagated instead of the 2(M x M) in Eq. (J|) and Eq. ([F]). 

III. STATE TO STATE APPROACH. 

The state to state approach || is based on the simultaneous optimization of N transitions 
between pairs of initial and final states. These pairs of states {<fu, ffi}, {I = 1,...,N) are 
related by the target unitary transformation, \(pfi) = O \(pu). The objective is formulated as 

N 

V (6;T;e) = Y,\(MT;e)\ ¥ > fl )\ 2 , (15) 
1=1 

where ipu(T\e) is the state at the target time that evolves with the driving field e(t) and 
the initial condition ipu(t = 0) = ipu. rj is a positive real number and its maximum value 
7] = N is reached when all the initial states ifn are driven by the field to the correct final 
states (ffi. The task of obtaining the optimal field is equivalent to the maximization of r\. 
The set of initial states (fin must be chosen carefully. In order to account for all the possible 
transitions the states ipn have to represent the relevant subspace. However, the choice of an 
orthonormal basis could produce undesired results. Let us denoted by {</?«} an orthonormal 
basis of the relevant subspace and by D an arbitrary unitary transformation diagonal in 
that basis. The product OD is also a unitary transformation. If eo and €od are fields that 
generate O and OD at time T respectively, it is found that the same optimization objective 
is obtained, 

J7j_(6;T;e )=7u(6;T;e OJJ ), (16) 

where _L denotes that rj was evaluated using an orthonormal basis. Then any algorithm 
based on t] and using an orthonormal basis would find the optimal field corresponding to 
any of the possible targets OD (O is a particular case when D is the identity operator). 
The reason is that r] is sensitive to the modulus of the projection of each pair {ipu, ip/i} but 
not to the relative phases between them. The phase problem can be overcome with a careful 
choice of the states (pn. Keeping the first N — 1 states of the basis and substituting the N 
state by YliLi 'Pu/VN, the maximum condition is achieved only when the field induces the 
target unitary transformation up to a possible global phase. 
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The optimization is formulated as the maximization of the functional 

K(il> ik ,i/> fk ,Ae) = y2\(MT;e + Ae)\v f i)\ 2 -\ [ \Ae(t)\ 2 dt 
, , Jo 

- N T 

J2(MT;e + Ae)\^ fl ) / (^| 



- 2 Re 



— + ^H(t; e + Ae) ) <-,,)(// 



(17) 



where ipu, ipfu e o an d Ae depend on time. The form in Ref. [[J is recuperated setting e = 
and interpreting Ae as the driving field. 

The variations with respect to ipu and ipif (I = 1, N) lead to the following 2N equations 

(18) 

with the condition \ipu(t = 0)) = \(fu), and 



— = ~^H(t; e + Ae)|^> 



-|^/,) = -^(t;eo + Ae)|V/i) 



(19) 



with the condition |V ; /fc(^ = T)) = \(ffk)- The variation respect to Ae leads to the equation 
for the correction to the field, 

' N 

Y,(^fi(T;e + Ae)\ Vfl )(^ fl {t;e + Ae)\fi\Mt;eo + Ae)) . (20) 



Ae ss (t) = -— Im 

A h 



1=1 

Eq. ( |TED and (ITS) are formally equivalent to Eq. (fTl|), propagating forward or backward 
on time the boundary conditions of the problem. Then the main difference between the two 
approaches is found in the expression for the correction to the field, Eq. @ and Eq. (p0|) . 
The connection between the two expressions is established by rewriting the states ip and ip 
as 

\<pa) -»■ 10. 
Wfi) - 6 10, 



|V«(*;e)> - U(t,0;e)|/), 
\Mt^)) - B\t,T;e)\l) 



(21) 



and Eq. 

Ae ss (t) 



^(/|U f (T, 0; e + Ae) 6|Z) (l\B{t, T; e + Ae) ft U(t, 0; e + Ae)|Z> 



(22) 
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This expression is formally equivalent to Eq. (|S|) except for the factor (1\\j\t, 0, e +Ae) 0\l). 
(The different factor 1/2 can be removed with a redefinition of A). Omitting this factor 
the two approaches would be completely equivalent: in that case the careful choice of the 
initial state set wouldn't be necessary and the two approaches would lead to the same set 
of equations. In general the equivalence is only formal due to the additional factor and 
the different set of states {\k)} and {\l}}- In the following we will chose \l) = \k) for 
I = 1, N - 1 and \l = N) = J2k=i \k)/VN. The condition for the optimal field (Ae = 0) 
is obtained using Eq. ( p0|) and (H), 



Im 



N 

^2(l\lj\T, 0; e opt ) 6\l) Jj\t, T; e opt ) £U(t, 0; e opt )\l) 



i=i 



0. 



(23) 



Eq. fl2"5| ) is a necessary but not sufficient condition for the optimal field, as is also Eq. ( |IUD . 
Denoting as before by e a driving field generating the target unitary transformation up to a 
global phase, the left-hand-side of (E3) gives, 



Im 



N 



^2(11^ v\t,T;e) p,U(t,T;e)6\l) 



0. 



(24) 



(The property Im[(\E f |/i| 1 I r }] = for Hermitian operator ft, was used). Contrary to Eq. (|i~2|). 
the phase <p doesn't appears in Eq. (p4]) due to the factor (1\\j\t, 0, e) O | Z) correcting each 
term in the sum. Using the same arguments leading to Eq. (|i~3|) , it can be shown that the 
left-hand-side of Eq. (p3|) is null for e = when the target unitary transformation is diagonal. 
The previous results imply that any field that generates the target unitary transformation 
up to a global phase fulfills the condition in Eq. (p3|) and in addition spurious solutions 
could be found, as in the first approach. 



IV. CONCLUSIONS 



In this study we have employed a modified formulation of OCT which uses the correction 
to the driving field Ae as central element. Ae = is a necessary but not sufficient condition 
for obtaining the objective. The relations derived from this condition allow a better analysis 
of the control equations, in particular, the similarities and differences between the two 
approaches. Optimal solutions of the evolution equation approach will be also optimal 
solutions of the state to state approach. In this sense both approaches are equivalents. The 
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difference between them is a term in the state to state approach which modified the phases. 
This term can cause a phase ambiguity in the target unitary transformation. A careful 
choice of the initial set of states can solve this problem. 
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